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Abstract 

The recent data for Bose-Einstein Correlations (BEC) of three-charged particles 
obtained by NA44 Collaboration have been analysed using theoretical formula with 
Coulomb wave functions. It has been recently proposed by Alt et al. It turns out 
that there are discrepancies between these data and the respective theoretical val- 
ues. To resolve this problem we seek a possibly modified theoretical formulation of 
this problem by introducing the degree of coherence for the exchange effect due to 
the BEC between two-identical bosons. As a result we obtain a modified formula- 
tion for the BEC of three-charged particles showing good agreement with the data. 
Moreover, we investigate physical connection between our modified formulation and 
the core-halo model proposed by Csorgo et al. Our study indicates that the interac- 
tion region estimated by the BEC of three-charged particles in the S + Pb collisions 
at 200 GeV/c per nucleon is equal to about 1.5 fm~1.8 fm. 

Key words: Bose-Einstein Correlation, three-charged particles, Coulomb wave 
functions, high energy heavy- ion collisions 



1 Introduction 



One of the most interesting subjects in high energy heavy-ion collisions is study 
of the higher order Bose-Einstein Correlation (BEC) effect [1-6] (known also 
as the HBT or the GGLP effect, or as the hadron interferometry [7-10]). From 
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data on BEC we can (in principle) infer the size of the interaction region and 
therefore estimate the energy densities reached in high energy collisions. Such 
work is a necessary task in the search for the quark-gluon plasma [11,12] - a 
new, hypothetical form of matter. 

To get more precise sizes of the interaction regions, we have to take into 
account the final state interactions among the charged particles [13,14]. A 
great advance in this direction for the BEC of the three-charged particles 
has been recently made by Alt et. al. [5]. They have derived a correction 
formula for the raw data introducing distribution functions of the charged 
particles. Their formulation is based on the plane wave functions and on the 
Coulomb wave functions, assuming that produced hadrons are already in the 
asymptotic region of the Coulomb interactions where the strong interaction 
already vanishes [15,16]. It amounts in the following correction factor K Cou i 
due to the Coulomb effect for identical three-charged particles: 3 

Kcoui = . (1) 



The denominator D p i ane is given by (p(xj) are distribution functions of charged 
particles): 



Dplane = ^ J rf 3 Xip(xi)rf 3 X2p(x 2 )rf 3 X3p(x 3 ) 



e i(ki-xi+k2-X2+k3-X3) _|_ e «(ki-X2+k 2 -xi+k3-X3) 
_|_ e i(ki-x2+k 2 -x 3 +k3-xi) _|_ e i(ki-xi+k 2 -x3+k3-x 2 ) 

2 

_|_ e i(ki-X3+k2-xi+k 3 -x 2 ) _|_ e i(ki-X3+k 2 -X2+k 3 -xi) ^) 



3 The correlation functions for two and three-identical particles are given as usual 
by 

^(2+ or 2-) p (k k ) ^(3+ or 3-) p (k k k ) 

and 



NBG P(ki)P(k 2 ) N B G P(k!)P(k 2 )P(k 3 ) ' 

where kj is the momentum of particle i, and P(ki, k2) and P(ki, k2, k 3 ) are two 
and three particles probability densities, respectively. The probability densities for 
two-identical particles case can be written as, 



// 



2 

^kik 2 ( x i> x 2) /5(xi)p(x 2 )d 3 xid 3 x 2 , 



where 



V>ki E k 2 ( X l> X 2 ) = -^=[^kik 2 (xi, X 2 ) + Vk 1 k 2 (x 2 , Xl)] 



p(xj) stand for the source functions of particle i. 
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The numerator Ncoui has the following form: 



N, 



Coul 



rf 3 Xip(x 1 )rf 3 X2p(x 2 )rf 3 X3p(x 3 ) 
V£k 2 ( X l, X 2)^£k 3 ( X 2> Xg)^!^. X l) 

+^k 1 k 2 ( x i ) x s)^k;k3( x 3 ? x 2)^ k ;k 1 ( x 2 ) x i) 
+^k;k 2 ( x 2, x i)^k 2 k 3 ( x i 5 x 3)^ k 7 3 k 1 ( x 3, x 2) 
+^ k 7 1 k 2 ( x 2, ^v^N, x i)^k 7 3 k 1 ( x i> x 2) 

+^£k 2 ( X 3, Xi)^k 2 k 3 ( X l> X 2)^£ i k 1 ( X 2, X 3 ) 

2 

+^k 1 k 2 ( x 3, x 2 )^£ 2k3 (x 2 , x i)^ k 7 3 k 1 ( x i ) x s) 



(3) 



Here V'k-k,-^' x j) are the Coulomb wave functions of the respective 2-body 
collision expressed as, 



^£ k .(Xi, x^rfl + ^je^^^'^Fh^, 1; i(fc T - - k^- • r^)],(4) 



with Tij — (Xj x jr')) ^ij — O^i K?)/^> '"ij — kij — Sid Tjij — 

ma/kij. F[a, b; x] and T(x) are the confluent hypergeometric function and 
the Gamma function, respectively. In order to use Eqs. (1), (2) and (3), one 
has to assume first some shapes and sizes for the source functions. In fact, this 
is the procedure already used in Ref. [17] by NA44 Collaboration: 

Corrected data = (raw data) x K spc x ^acceptance x K Cou i , 

where K spc and K acceptance denote the effect of multiparticle production in the 
single particle spectra and the acceptance effect in the experiment. 

In this paper, we would like to adopt a different point of view for Eq. (3). As is 
seen in Ref. [14], the BEC of identical two-charged pions can also be analysed 
by the Coulomb wave functions. It is therefore reasonable to expect that the 
numerator Ncoui is the main theoretical ingredient in analysis of the BEC of 
three-charged particles. We argue therefore that 

^(3+ or 3-) /N BG ^ Cy< ? (5) 



where we have introduced the normalization factor C, which corresponds to 
the asymptotic value of the BEC. Using Eq. (5) we can now (with the help of 
the CERN-MINUIT program) analyse data of Ref. [17] using Gaussian source 
distributions of radii R, p(x) = j^Wj^ ex P — W 1 4 ■ 

4 It should be remembered that NA44 Collaboration data are for the variable 



Q§(4D) = (h - k 2 ) 2 + (k 2 - k 3 ) 2 + (fc 3 - fci) 
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In the next paragraph, we analyse the data of NA44 Collaboration [17] by Eq. 
(5). In the third paragraph we shall derive a modefied theoretical formula for 
3-particle BEC introducing the degree of coherence parameter into Eq. (5). 
This formula will be then used in the 4th paragraph for the re-analyses of the 
experimental data [17]. In the 5th paragraph, we investigate whether or not 
there is physical connection between our study and the core-halo model [18]. 
Concluding remarks are given in the final paragraph. 



2 Application of Eq. (5) to the data by NA44 Collaboration 



Here we analyse the data by Eq. (5). As can be seen in Fig. 1 and Table 1, 
there are some discrepancies between the data points and theoretical values 
calculated by means of Eq. (5). The minimum x 2 value, 17.6 in Table 1, seems 
to be big, as the number of the data points are considered. Thus we would 
like to know why this equation cannot explain the data [17]. We know that 
there are several possible reasons due to effects of the partial coherence, the 
contamination [19] and the long-lived resonances [18]. At present, we consider 
the effect of the possible partial coherent of produced pions. In fact, authors 
of Ref. [17] have used not the equivalence of Eq. (5) but the following formula 
instead (cf., Ref. [9,20]): 

!*£ = C(1 + *-**). (6) 



It contains one more parameter, A 3 , which can be regarded as a kind of effective 
degree of coherence and which, in our opinion, should therefore occur also 
somehow in Eq. (5). 



where ki are four-momentum of charged particles. Q3 = yQ%- However, in our 
calculations we assume that q^ ^ = (k j — k .,) 2 ~ an d use, instead, 

Q 2 3 (3D) = - k 2 ) 2 + (k 2 - k 3 ) 2 + (k 3 - kx) 2 . 
There is an approximation between them 

Q 3 (AD) * Q 3 (3D) - 5>o,ii/2Q3(3£>) . 

Qs(AD) depends on sum of squared energy differences. 
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Fig. 1. Analysis of 3ir + BEC in S + Pb collision [17]. (a) and (b) are results of Eqs. 
(5) and (6), respectively. The error bars are sum of statistical and systematic errors. 

Table 1 

Estimated values for the data [17] by Eqs. (5) and (6) using CERN-MINUIT pro- 
gram. 



Formulas 


C 


R [fm] 


A 3 


X 2 /N dof 


Eq. (5) 


0.941±0.026 


2.47±0.14 




17.6/16 


Eq. (6) 


0.986±0.028 


2.36±0.26 


1.37±0.19 


7.8/15 



3 Diagram Decomposition of Eq. (5) 



First of all, we have to find a possible way for the introduction of the degree of 
coherence parameter A into Eq. (5). Let us therefore examine the plane wave 
(PW) approximations of the Coulomb wave functions: 



^( 1 )=^k 1 k 2 ( X l> X 2)^k 7 2 k 3 ( X 2> XS^kjfXS, Xi) 

p W e «ki2-ri2 ik23-r23gjk3i-r3i _ (3/2)j(ki -xi+k2-X2+k 3 -X3) 

^(2)=<k 2 ( X l, X 3Xk 3 ( X 3, X 2 )V£ kl (x 2 , Xl ) 

PW , e iki2Ti3 e ik23T32 e jk 3 iT2i _ e (3/2)i(ki-xi+k 2 -X3+k 3 -X2) 

^(3)=^£k 2 ( X 2, x l)^k 2 k 3 ( X l> X 3)^£ i k 1 ( x 3, x 2) 

PW , giki2T2i e ik2 3 Ti3 e ik3iT32 _ e (3/2)j(ki-x 2 +k2-xi+k 3 -X3) 

^(4)=<k2( x 2> x 3 )< k3 (x 3 , x^g^Xi, x 2 ) 

PW ) giki2T23gik23T3i e jk3iTi2 _ g (3/2)i(ki ■x 2 +k2-x 3 +k 3 -xi) 

^(5)=^k 1 k 2 ( x 3, x l)^k 7 2 k 3 ( x l ) x 2)V' k ;k 1 ( x 2 ) x 3) 

p W e iki2-r3i e ik2 3 -ri2gik3i-r23 _ e (3/2)i(ki-x 3 +k2-xi-)-k 3 -X2) 

^(6)=<k 2 ( x 3> x 2 )^ k3 (x 2 , x^g^Xi, x 3 ) 
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p W giki2T32 e ik23T2i e ik 3 iTi3 _ e (3/2)i(ki-x 3 +k2-X2+k 3 -xi) 



(7f) 
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Fig. 2. Diagram reflecting three-charged particles Bose-Einstein Correlation (BEC) 
and Coulombic potential (V^). X means the exchange effect of BEC. 

Notice that, except for the factor 3/2, exponential functions are the same 
expressions as those present in the integrand of Eq. (2). This difference is 
attributed to the fact that Coulomb wave function used here describes two- 
charged particles collisions, therefore factor 3/2 appears because there are 
relevant two-particle three combinations among three-charged particles. 



Combining Eqs. (7) and Figs. 2, we obtain the following three sets of equations: 



F 1 = ^#)A1^1, (8a) 

F 2 = ^[A(1)A*(2) + A(1)A*(3) + A{1)A*{6) + A(2)A*(A) + A{2)A*(h) 
+A(3)A*(4) + A{3)A*{5) + A(4)A*(6) + A(5)A*(6) + c. c] 

BEC between two-charged particles (8b) 
(See Figs. 2 (b)~(d)) , 

F 3 = i[A(l)A*(4) + A(1)A*(5) + A(2)A*(3) + A(2)A*(6) + A{3)A*{6) 
+A(A)A*(5) + c. c] 

— > BEC among three-charged particles (8c) 
(See Figs. 2 (e) and (f)) . 
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Combining now Eqs. (8) and the concept of partial coherence for the BEC 
[17,20], we can introduce a coherence parameter v^A for the single mark (x) 
in Fig. 2. The A = 1 corresponds to the totally chaotic source, which is the 
assumption behind Eq. (5). Taking into account the strength of the degree 
of coherence A between two-identical bosons and A 3//2 among three-identical 
bosons in Figs. 2, we can finally express the BEC for three identical charged 
particles as: 

at(3+ or 3-) , 

NBG =C j rf 3 x lP (x 1 )rf 3 x 2 p(x 2 )rf 3 X3p(x 3 )[F 1 + XF 2 + A 3 / 2 F 3 ] . (9) 

Equation (9) is the modified theoretical formula we were looking for. It differs 
from Eq. (5) originally proposed by Alt et. al. in [5] by the presence of the 
degree of coherence A and in the limit of rjij — > it becomes 

Eq . (9) C (l + 3Ae"f + 2A 3 /Vf ^) , (10) 

which is the extended formula proposed some time ago by Deutschmann et al. 
[20]. 

It should be noticed that Eq. (9) can be applied to data corrected only by 
the Gamow factor G{rji2)G{rj2z)G{rj^i) in an ideal case [21], because Eq. (9) is 
described by the Coulomb wave functions including the Gamow factors (see 
Ref. [22]) 5 . 



4 Reanalyses of NA44 Collaboration data by means of Eq. (9) 

At present we have no data corrected only by the Gamow factors, therefore 
we apply Eq. (9) to the analysis of NA44 Collaboration data [17] using the 
CERN-MINUIT program. Our results are shown in Fig. 3 and Table 2. Com- 
paring them with those of Table 1, it can be said that the x 2 -value becomes 
smaller, i.e., the agreement is now better. The range of interaction becomes 
also smaller. For the sake of reference we present in Table 2 also results ob- 
tained by using Eq. (10) . 



5 In other words, ideal data sets for Eq. (9) are of the form of 

Corrected data = (raw data) x -Kgpc x -^acceptance amow j 

where K Gamow = 1/ {G(f]i 2 )G{r]22,)G{r] Z1 )). 
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Fig. 3. Reanalyses of 3vr + BEC in S + Pb collision [17]. (a) is result of Eq. (9). (b) 
is that of Eq. (10). 



Table 2 

Reanalyses of 3tt+ BEC in S + Pb collision [17] by Eqs. (9) and (10). 



Formula 


C 


R [fm] 


A X 2 /N dof 


Eq. (9) 


0.917±0.032 


1.53±0.20 


0.55±0.07 6.7/15 


Eq. (10) 


0.984±0.029 


2.60±0.28 


0.33±0.04 7.7/15 



From results of Table 2, we see that R— 1.53 fm by Eq. (9) is small. The reason 
is attributed to the fraction of partially coherent effect (A). From comparisons 
between results by Eqs. (9) and (10), it can be seen the interaction region 
becomes smaller, and the degree of coherence does conversely bigger due to 
the Coulomb interaction. 
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5 Possible interpretation of A by core-halo model 



The core-halo model has been proposed by Csorgo et al. [18]. We study whether 
or not there is physical connection between our previous formulation and theirs 
[18]. To apply their model to Eq. (8), first of all we use the following source 
functions 

P(Xi) = Pc(Xi) + phalo(Xi) , (11) 



where p c (xi) = 




Fig. 4. Physical picture of core-halo model. 



In concrete analyses, the radius of the halo becomes to be infinite, due to 
the effect of long-lived resonances (Rh — > oo). The exchange function of the 
halo part due to the BEC, E ha i Q = exp(—Q 2 Rl) — > 0. Then we can interpret 
A in Eq. (9) by their parameter, i. e., fraction of multiplicity from the core 
part, f c = (n core ) / (n to t) , as A = f%. In this case, it should be stressed that 
their parameter p c = (n co ) / (n core ) , the fraction of the coherently produced 
multiplicity from the core part is zero. 

Moreover, we can assume the laser optical approach for the core part [18], 
introducing a parameter p = (n c / iao )/(n core ) [1]. Notice that p = 1 — p c [18]. 
For the cross mark (x) in Fig. 2, we assume two components of coherently 
(1 — p) and chaotically produced particles (p). Then we obtain the following 
expression 

J rf 3 x lPc (x 1 )rf 3 x 2 p c (x 2 )rf 3 X3p c (x 3 )(F 1 + f c p 2 F 2 + / c 3 p 3 F 3 ) 
+ / d 3 xip c ( Xl ) • 5 3 (xi)d 3 x 2 p c (x 2 )d 3 X3p c (x 3 ) 



yy(3+ or 3-) 
N BG 
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• (2/ c 2 p(l - p)F 2 + 3/y (1 - p)F 3 ) } . (12) 
The effective degree of coherence, A3, the intercept at smallest Q 3 is given as 
A* = / 2 (p 2 + 2p(l - p)) + f c (p 3 + 3p\l - p)) . (13) 



In Eq. (12), as p = 1, we obtain Eq. (9) with A = / 2 . On the contrary, as 
f c — 1, we obtain an expression of laser optical approach [1]. By making use 
of Eq. (12) and an expression for the two-charged particles of the BEC, we 
can analyse the data by NA44 Collaboration. 

For the two-charged particles of the BEC by means of the core-halo model 
with laser optical approach, we have the following equation, 

7^(2+ or 2-) r , 



Results of our analyses are shown in Table 3 and Figure 5. As is seen in 
them, the interaction ranges of i?(core part) are estimated in the ranges of 
1.5 fm < i?(core) < 1.8 fm, and 4.7 fm < i?(core) < 5.4 fm, for the BEC of 
3n and 3n — > 2ir, respectively. 

The common region between results of the BEC of three-pion and two-pion is 
roughly described by 0.6<p<1.0 and f c ~ 0.7. 



6 Concluding remarks 

We have derived the theoretical formula for the BEC of three-charged identical 
particle using both the Coulomb wave functions and the notion of the degree 
of coherence and compared it with the experimental data 6 . Historically the 

6 For the numerical calculations of Eqs. (5) and (9) (in order to save the CPU-time), 
we have first calculated 10 3 k values of the Coulomb wave functions, which were then 




•■0kik 2 ( x 2> x i)j- The effective degree of coherence A?; is given as 



A; = / c V + 2p(l-p)). 



(15) 
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Table 3 

Typical results from analyses of data of NA44 Collaboration by Eqs. (12)~(15). 



3vr BEC 


V 


1.0 


0.8 


0.6 


fc 


0.743 ± 0.050 


0.779 ± 0.055 


0.856 ± 0.061 


A3 


0.964 


1.007 


1.022 


R (fm) 


1.53 ± 0.21 


1.72 ± 0.25 


1.87 ± 0.27 


X 2 /N dof 


6.7/15 


6.6/15 


6.6/15 




3vr -f 


2tt BEC 




P 


1.0 


0.8 


0.6 


fc 


0.633 ± 0.028 


0.647 ± 0.029 


0.688 ± 0.032 


A2 


0.400 


0.402 


0.398 


R (fm) 


4.69 ± 0.45 


5.34 ± 0.54 


5.85 ±0.59 


X 2 /N dof 


14.6/17 


14.9/17 


14.9/17 




Fig. 5. Sets of f c and p estimated in analyses of BEC of two-charged and 
three-charged pions using the CERN-MINUIT program. Widths of / c 's stand for 
error bar of ±2<r. 

degree of coherence in the BEC of the two-identical bosons has been intro- 
duced by experimentalist [20], and theoretical works in this direction have 
been performed in Ref . [9] . 

Our first analyses suggest that the degree of coherence A is a necessary in- 
gredient also for the BEC of three-charged particles, in the same way as it 
was for the BEC for two-charged identical particles. This fact means that the 



used together with some interpolation procedure during the concrete calculations. 
In this way we could make use of the CERN-MINUIT program in our analyses. 
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source producing finally observed particles is not purely chaotic. However the 
investigation in 5th paragraph elucidates that our modified formulation can 
be interpreted by the core-halo model. Our degree of coherence A is equal to 
fl [introduced in Ref. [18]], provided that particles are chaotically produced; 
A = / c 2 (= ((n core ) I (n tot )) 2 ) . Moreover, if we can assume the laser-optical ap- 
proach for the cross mark ( X ) in Fig. 2, we obtain Eqs. (12) ~ (15). By making 
use of them, we obtain Table 3 and Fig. 5. There is common region among 
results from analyses of the BEC of 3n and 3n — > 2n processes, 0.6^p^l.O 
and f c ~ 0.7. However, the magnitude of the interaction regions estimated by 
the BEC of 3n and 3n — > 2n process are different. This problem should be 
considered in the future. 

It should be noticed that also recent data on the BEC of 3n~ reported by 
OPAL Collaboration [23] suggest the necessity of introduction of some degree 
of coherence or f c in Ref. [18] 7 . 
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Addendum 



Because we had no sufficient information on the raw data and corrected data in 
2001, we utilized the corrected data with -Kcouiomb in our analyses. Of course we 
mentioned that our Eqs. (9), (12), and (14) are available for the corrected data 
including -Kspc (single particle correction (SPC)) and -^acceptance, explicitly. 
Referring to Refs. [25,26], we have examined the methods of correction used 
by NA44 Collaboration [17]. In this addendum, thus we can present the data 
including K SPS and -^acceptance and analyze them by means of Eqs. (9), (12) 
and (14). Before concrete analyses, we categorize raw and corrected data in 
Table 4. 
Table 4 

Category of data. Notice that R = input is necessary for -Kcouiomb, where R = 5 fm 
is used in Refs. [17,25,26]. 

1) Raw data raw data 

2) Quasi-corrected data (Q-CD) (raw data) x -fTsPC x -^acceptance 

3) Corrected data with Gamow (raw data) x -fTsPC x ^acceptance x -Kcamow 

4) Corrected data with Coulomb (raw data) x ifsPC x ^acceptance x -^Coulomb 



Using the data in Ref. [17] with -ft'couiomb or -ft'camow, we can obtain the quasi- 
corrected data (Q-CD) shown in Fig. 6, which are available for our purpose. 
Using Eqs. (5) and (9), we obtain the results shown in Table 5. As seen in 
Fig. 6, the our formulation including the degree of coherence A, i.e., Eq. (9), 
seems to be available for analyses of the charged ir BEC. 




03 ( MeV > 

Fig. 6. Analyses of quasi-corrected data (Q-CD). (a) and (b) are results of Eqs. (5) 
and (9), respectively. . 

To examine availability of the core-halo model described by the Coulomb wave 
functions, i.e., Eqs. (12) and (14), we analyze the quasi-corrected data (Q-CD). 
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Table 5 

Estimated values for the data [17] by Eqs. (5) and (9) using CERN-MINUIT pro- 
gram. 



Formulas 


C 


R [fm] 


A 


X 2 /N do f 


Eq. (5) 


0.94±0.03 


5.83±0.40 




30/16 


Eq. (9) 


0.91±0.03 


2.89±0.39 


0.45±0.05 


6.7/15 



Our results are shown in Fig. 7 and Table 6. Coincidences between 3tt BEC 
and 3ir — > 2tc BEC in Fig. 7 ( f c vs. p) seem to be more enlarged than those 
of Fig. 5. 

In this addendum, we have analyzed the quasi-corrected data (Q-CD) named 
in Table 4, utilizing Eqs. (5), (6), (12) and (14). It can be said that our 
formulation works well for the description of the charged 3n BEC. 

1 

0.8 



0.6 
0.4 
0.2 


0.2 0.4 0.6 0.8 1 

P 

Fig. 7. Sets of f c and p estimated in analyses of BEC of two-charged and 
three-charged pions using the CERN-MINUIT program. Widths of / c 's stand for 
error bar of ±2a. 

As seen in the explanation above, Fig. 1(a) and the upper line (Eq. (5)) of 
Table 1 , and Fig. 3(a) and the upper line (Eq. (9)) of Table 2 should be 
replaced by Fig. 6 and Table 5. Moreover, Fig. 5 and Table 4 should be replaced 
by Fig. 7 and Table 6. 

One of authors (M. B.) would like to thank H. B0ggild, T. Csorgo and B. 
Lorstad for useful conversations. 
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Table 6 

Typical results from analyses of data of NA44 Collaboration by Eqs. (12)~(15). 
Rpinne = (3/2)i?, where i? p i ane is corresponding to a parameter in the plane wave 
formulation. 



p 


1.0 


0.8 


0.6 




3vr BEC 




fc 


0.67±0.04 


0.70±0.04 


0.76±0.05 


AS 


0.75 


0.77 


0.78 


R (fm) 


2.89±0.39 


3.22±0.46 


3.47±0.51 


-Rplane (f m ) 


4.33±0.58 


4.83±0.69 


5.21±0.76 


X 2 /N dof 


6.7/15 


6.8/15 


6.8/15 




3^ 2tt BEC 




fc 


0.67±0.03 


0.69±0.03 


0.73±0.03 




0.45 


0.46 


0.45 


R (fm) 


4.42±0.38 


5.01±0.45 


5.49±0.50 


X 2 /N dof 


14.5/17 


14.7/17 


14.7/17 
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